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In quantum chromodynamics (QCD), the binding energy of an infinitely heavy quark-antiquark pair in a 
color singlet state can be calculated as a function of the distance. We investigate this static potential of QCD 
perturbatively and calculate the full two-loop coefficient, correcting an earlier result. Beyond this order, the 
perturbative expansion breaks down. 



1. Introduction 

The static potential of quantum chromody- 
namics (QCD) is subject to theoretical inves- 
tigations since more than twenty years. Be- 
ing the non-abelian analogue of the well-known 
Coulomb potential of quantum electrodynamics 
(QED), this interaction energy of an infinitely 
heavy quark-antiquark pair is a fundamental con- 
cept which is expected to play a key role in the 
understanding of quark confinement. Moreover, 
the static potential is a major ingredient in the 
description of non-relativistically bound systems 
like quarkonia, and it is of importance in many 
other areas, such as quark mass definitions and 
quark production at threshold. 

It is expected that the static potential con- 
sists of two terms: a Coulomb-like term at short 
distances, which is calculable with perturbative 
methods, and a long-distance term responsible 
for confinement. Even though a perturbative 
analysis is not suited to give the full potential, 
such a calculation proves very useful. The short- 
distance part of the potential can be utilized as a 
refined starting point for the construction of po- 
tential models (which have been rather successful 
in the past for the description of quarkonia), or 
it could describe very heavy systems (like it) to 
good accuracy. Furthermore, it can be compared 
to the results of numerical calculations in lattice 
gauge theory. It is natural to define the QCD 
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coupling constant with help of the potential as 
V(r) = _| ^ tne so _ ca n ec i V-scheme, us- 

ing a physical quantity in contrast to the usual 
coupling definition in the MS scheme. In lattice 
calculations a v is regarded as as the 'better' ex- 
pansion parameter [ For these reasons, and 
to get a more precise determination of a—^ from 
the lattice, the relation between the two couplings 
has to be known. 

A first determination of the static poten- 
tial in (massless) QCD has been performed by 
L. Susskind in the context of a lecture about lat- 
tice gauge theory [ 0. In order to demonstrate 
asymptotic freedom in Yang-Mills theory, he cal- 
culated the one-loop pole terms using a Wilson- 
loop formula for the potential, and re-derived 
the first coefficient of the renormalization group 
Beta function. This work was extended by other 
groups quite soon, who then added fcrmionic con- 
tributions [ H and two-loop pole terms [ || to the 
potential, as well as examined the structure of 
higher-order corrections qualitatively [ |[. Re- 
cently, the perturbative static potential has re- 
ceived new interest, in particular due to its appli- 
cation in top-quark production at threshold [ || , 
a process that comes within experimental reach in 
the near future. A complete two-loop calculation 
for the static potential was performed in ref. [ Q . 
Such an important result clearly needs confirma- 
tion. This is one of the motivations of our work 
on the two-loop potential [ |l |||. More recently, 
the effect of fcrmion masses was considered on the 
two-loop level in ref. [ [h]] . 
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2. Definition and Expansion 

The static potential is denned in a manifestly 
gauge invariant way via the vacuum expectation 
value of a Wilson loop [ p|, pd||, 

V(r) = - hm iln(Wr) , (1) 

1 — *oo 1 

W r = trPexp (ig j> dx^A^j . (2) 

Here, T is taken as a rectangular loop with time 
extension T and spatial extension r. The gauge 
fields are path-ordered along the loop, while 
the color trace is normalized according to tr(..) = 
tr(..)/trl . 

In a perturbative analysis it can be shown that, 
at least to the order needed here, all contributions 
to eq. (|l]) containing connections to the spatial 
components of the gauge fields Ai(r,±T/2) van- 
ish in the limit of large time extension T. Hence, 
the definition can be reduced to 

Vrt = ~ r Hm ^ In (trT exp (jjt A l) ) > ( 3 ) 

where T means time ordering and the static 
sources separated by the distance r = \r — r'\ 
are given by 

J°(x)=ig5 fi0 T a [S(x-r)-S(x-r')} , (4) 

where T a are the generators in the fundamental 
representation. In the case of QCD the gauge 
group is SU(3). The results will be presented 
for an arbitrary compact semi-simple Lie group 
with structure constants defined by the Lie al- 
gebra [T Q ,T b ] = if abc T c . The Casimir oper- 
ators of the fundamental and adjoint represen- 
tation are T a T a = C F and f acd f bcd = C A S ab . 
tr(T a T b ) = T F 5 ab is the trace normalization, 
while n/ denotes the number of massless quarks. 

Expanding the expression in eq. (^) perturba- 
tively, one encounters in addition to the usual 
Feynman rules the source-gluon vertex igS^oT 11 , 
with an additional minus sign for the antisource. 
Furthermore, the time-ordering prescription gen- 
erates step functions, which can be viewed as 
source propagators, analogous to the heavy- 
quark effective theory (HQET). 
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Figure 1. Classes of two-loop diagrams con- 
tributing to the static potential. Double, wiggly, 
dotted and solid lines denote source, gluon, ghost 
and (light) fermion propagators, respectively. A 
blob on a gluon line stands for one-loop self- 
energy corrections. 



Concerning the generation of the complete set 
of Feynman diagrams contributing to the two- 
loop static potential, there are some subtleties 
connected with the logarithm in the definition (||). 
All this is explained in detail in [ 0, , so we 
only list the relevant exchange diagrams here (see 
fig-0). 

3. Method 

Turning to the purely technical side of the 
work documented here, only part of the Feyn- 
man diagrams to be considered, namely the pure 
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self-energy contributions to the static potential, 
are amenable to standard calculation methods 
[ ||, H II IJ. For the others, essentially 
being two-point functions also, but involving 
non-covariant propagators, the standard meth- 
ods need to be generalized. We have developed 
a general strategy to deal with these expressions, 
which is based on a purely algebraic reduction to 
a minimal set of integrals [ || . 

The method employed can be briefly summa- 
rized as follows: 

• Working in momentum space, all dimen- 
sionally regulated (tensor-) integrals are re- 
duced to pure propagator integrals by a gen- 
eralization of the method of T-operators 
[ |l4|j . The resulting expressions are then 
mapped to a minimal set of five scalar in- 
tegrals by means of recurrence relations, 
again generalizing [ [TJ|] as well as [ |l^| to the 
case including static (non-covariant) prop- 
agators. These two steps have been imple- 
mented into a FORM package. Thus, we 
constructed our method to be complemen- 
tary to the calculation in ref. [ assuring 
a truly independent check of the result pre- 
sented therein. At this stage, one obtains 
analytic coefficient functions (depending on 
the generic space-time dimension D as well 
as on the color factors and the bare cou- 
pling), multiplying each of the basic inte- 
grals. 

• The basic scalar integrals are then solved 
analytically. Expanding the result around 
D — 4 (which is done in both MAPLE and 
Mathematica considering the complexity of 
the expressions), renormalizing and Fourier 
transforming back to coordinate space, one 
obtains the final results presented below. 

Important checks of the calculation include 

• gauge independence of appropriate classes 
of diagrams, 

• confirmation of cancelation of infrared di- 
vergences, 

• correct renormalization properties. 



4. Result and Discussion 

We obtain, as our final result for the two-loop 
static potential in coordinate space, 

a v (l/r 2 ) 



V(r) = -C F ■ 



with 



(5) 



a :_ 



a v {l/r 2 ) = i + _^( ai+/ 5 L) 
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+ TeS I ° 2 + ™ + — ) + (ft + 2 ^i)l 



(6) 



where L = 2ln(ure 7 ) and a = a (u 2 ). The 

' MS ^ MS Kr ^ ' 

first two terms of the beta function are given 
by /3b = ^ C A - | T F n f and ft = f G\ - 
AC F T F nf - ^ CAT F nf. The one- and two-loop 
constants read 
31 
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respectively. As it has to be, the coefficients prove 
to be gauge independent. Comparing our two- 
loop result for 02 with [ [?j , we find a discrepancy 
of 2tt 2 in the pure Yang-Mills term (oc C\). This 
amounts to a 30% decrease of a-i for the case of 
rif = 0, and a 50% decrease for nj — 5 (for 
SU(3)), which is the case needed for ti thresh- 
old investigations. This difference can be traced 
back to a specific set of diagrams, and it is clari- 
fied with the author of ref. [ 0] 

There are now numerous concepts for a 'renor- 
malization group improvement', i.e. for an 'opti- 
mal' choice of the MS scale parameter /2 in order 
to reduce large logarithmic corrections. Exam- 
ples include the 'natural choice' /2 = e~ 7 /r and 
the choice p, — exp(— 7 — ai/2/3o)/r, which elim- 
inates the one-loop coefficient completely. Due 

2 Wc thank M. Peter for checking this result. 
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to this freedom, it is not very illuminating to 
present plots of the coordinate space potential. 
A general feature is that, at increasing distance, 
the large two-loop coefficient begins to dominate 
quite soon, even causing the potential to decrease 
again above some r cr it, to signal that the pertur- 
bative approach can be followed up to this critical 
distance at most. For a discussion of a variety of 
scale choices we refer to ref. [ Q . The smaller coef- 
ficient a 2 found in our calculation does not change 
the plots presented there qualitatively. The rea- 
son is that the term ^-(3q, which shows up as 
a result of the Fourier transform, is of the same 
order as 02. 

Concerning the convergence of the series (^), 
let us give some numbers. For SU(3), Tp = \ 
and the 'natural' scale choice, we have 

a v = a m { 1+^(2.6-0.3^) 

+ (53.4 - 7.2n f + 0.2n 2 f ) + .. U9) 

In ref. [ , the first number in the two- loop term 
was 64.5. Apparently the convergence properties 
do still not look very promising. 

A comparison with four-dimensional quenched 
lattice results is given in ref. [ [l(|. There, it 
is concluded that the perturbative potential al- 
ready fails to describe the slope of the lattice po- 
tential at the smallest distances that are numeri- 
cally tractable, hence invalidating the possibility 
to match the two potentials at small distances. 

Summarizing, we have re-calculated the two- 
loop static potential by a method complemen- 
tary to the approach in [ f?J. We have devel- 
oped an algorithm which enables us to work in 
general covariant gauges throughout. Confirm- 
ing the fermionic contributions to the two-loop 
coefficient 02, we find a substantial deviation in 
the pure gluonic part of ai. The source of the 
discrepancy could be identified. The bad con- 
vergence of the perturbative series does not im- 
prove considerably taking into account the new 
value of a.2- Hence, the use of a physical cou- 
pling, defined by the potential, as expansion pa- 
rameter seems to be disfavored. Further studies 
are clearly needed to clarify the role of higher- 



order corrections. There has always been discus- 
sion about whether the perturbative Wilson-loop 
formula is a good definition of the static potential, 
which can be questioned due to possible infrared 
divergences in higher orders [ [5|]. A redefinition 
was proposed recently [ [Tt]] , which becomes effec- 
tive at the three-loop level. 
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